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Abstract 

A theoretical study of toroidal membranes with various degrees of intrinsic 
orientational order is presented at mean-field level. The study uses a simple 
Ginzburg-Landau style free energy functional, which gives rise to a rich vari- 
ety of physics and reveals some unusual ordered states. The system is found 
to exhibit many different phases with continuous and first order phase tran- 
sitions, and phenomena including spontaneous symmetry breaking, ground 
states with nodes and the formation of vortex-antivortex quartets. Transitions 
between toroidal phases with different configurations of the order parameter 
and different aspect ratios are plotted as functions of the thermodynamic pa- 
rameters. Regions of the phase diagrams in which spherical vesicles form are 
also shown. 
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I. INTRODUCTION 



The bilayer fluid membranes, which can form spontaneously when molecules with hy- 
drophobic and hydrophilic parts are introduced to water, are a popular topic of research. In 
theoretical studies of such membranes, it is usual to work on length scales much larger than 
the molecular size, so that they can be considered as continuous surfaces; two-dimensional 
curved spaces embedded in Euclidean three-space. Features commonly incorporated into 
mathematical models of membranes include bending rigidity k (disafiinity for extrinsic cur- 
vature), Gaussian curvature modulus (disaffinity for intrinsic curvature), spontaneous 
curvature due to membrane asymmetries, and for closed surfaces, constraints of constant 
volume, constant area and constant difference in area between the inner and outer layers. 
In recent years, interest has grown in membranes whose molecules have orient at ional order 
within the surface. For instance, a smectic-A liquid crystal membrane, whose molecules have 
aliphatic tails which point in an average direction normal to the local tangent plane of the 
membrane, can undergo a continuous phase transition to smectic-C phase, in which the tails 
tilt. The tilt is described by a two-component vector order parameter within the surface, 
given by the local thermal average of vectors parallel to the tails, projected onto to local 
tangent plane. If such a membrane has intrinsic (Gaussian) curvature, its orientational order 
is frustrated, since a parallel vector field cannot exist in a curved space. Ordering is further 
frustrated on a closed surface of spherical topology (genus zero), as an intrinsic vector field 
must have topological defects whose indices sum to the Euler number of the surface; in this 
case, two. ie. there must be at least two vortices in the order parameter field on a sphere. 
The order parameter field vanishes at these defects to avoid infinite gradients. The defects 
are energetically costly but topologically unavoidable. The energetics of a smectic-C, genus 
zero membrane were studied by MacKintosh and Lubensky in |l|]. There are other types of 
orientational order. Nematic order is similar to vector order except that its order parameter 
is locally invariant under rotations through 180°. Hence nematics can form defects of index 
i. Hexatic order is locally invariant under rotations through 60°, as this is the molecular 
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'bond' angle. The notion may be generalized to fluid phases with n-atic orientational order 
being locally invariant under rotations through ^ and forming defects of index ^- If a vector 
('monatic') is represented by an arrow then we may picture an ra-atic as an ra-headed arrow 
or an n-spoked wheel. For a theoretical study of n-atic fluid membranes of genus zero, see 

0- 

In this paper, theoretical findings on n-atic fluid membranes of toroidal topology are 
presented. As the torus has genus one and Euler number zero, it is possible to arrange an 
intrinsic vector field on a torus without topological defects, {ie. furry tori can be brushed 
smooth.) Therefore it appears plausible that n-atic membranes with large order parameter 
coupling compared with bending rigidity may energetically favour a toroidal topology to 
a spherical one. The object of embarking on this research was to confirm or refute this 
hypothesis. As will be demonstrated, it turns out to be true, although the arrangements of 
ordering in the toroidal surface are far more diverse than expected. 

Throughout this paper, a number of approximations are used. One approximation is the 
use of mean field theory. That is, no thermal fluctuations are considered. So the configura- 
tion adopted by the system is assumed to be that which minimizes its free energy. Also, the 
lowest Landau level approximation is employed in constructing the order parameter field, 



as discussed in section The validity of these approximations is investigated in section 
^ Both of these approximations give rise to apparent phase transitions between separate, 
clearly defined states of the system, some of which may in fact be blurred by thermal fluc- 
tuations and higher Landau levels, to the point where they disappear. It should therefore 
be understood that, where reference is made to phase transitions, it is in the context of 
mean field theory and the lowest Landau level approximation. Before studying the thermal 
physics, it will be necessary to introduce some tensor calculus. This is needed to describe 
the behaviour of an n-atic field living within a two-dimensional curved space, as well as 
the intrinsic and extrinsic elastic properties of the space {ie. the membrane) itself. In the 
absence of intrinsic ordering {ie. under the influence of the membrane's elasticity alone), one 
might expect the equilibrium shape of a toroidal vesicle to be axisymmetric. However, such 



a vesicle is only neutrally stable with respect to one particular mode of non-axisymmetric 
deformation. This Goldstone mode is discussed in section [V|, in which a further approxima- 
tion is introduced, limiting the range of vesicle shapes to be explored. In the sections that 
follow, the ground-state configurations and energies of the order parameter field are calcu- 
lated for each of these shapes, subject to the lowest Landau level approximation. Although 
some of these shapes are never energetically favoured, the behaviour of the field living on 
them is of interest, if only academic. Finally, in section |X|, the total free energy (consisting 
of order-parameter field and membrane elasticity parts) is minimized with respect to vesicle 
shape, to find the equilibrium configuration of the system as a function of temperature and 
the other thermodynamic parameters. The result is a set of phase diagrams for various 
values of n, which contain phases of vesicles with different shapes and ordering. The shapes 
should be experimentally observable, although the intrinsic ordering is most likely not. 



II. VALUES OF N 

n is the order of rotational symmetry. As stated earlier, the fluid membrane is locally 
invariant under rotations through and the n-atic field may be represented by a set of 
n-spoked wheels tessellating. Clearly n is integer. Furthermore, it must be one of the set of 
numbers for which regular n-gons tessellate. That is, n G {1, 2, 3, 4, 6}. These are the values 
of n which Park et. al. studied in 0]. However, consider a triangular {n = 3) lattice with 
long-range orientational, but no positional order; the n = 3 analogue of a hexatic lattice. 
Bonds between neighbouring molecules lie at angles which are integer multiples of |. Its 
symmetry properties in the absence of positional order are therefore indistinguishable from 
its dual lattice, the hexatic, and it will form defects of index |, in preference to the more 
energetically expensive index | defects. To convince oneself of this, it might be helpful 
to attempt to draw a lattice which has 120° rotational symmetry, but does not have 60° 
rotational symmetry. This is only possible if one keeps strict account of the positions of 
the lattice sites and polygons {ie. in the absence of dislocations). The following research is 



4 



therefore conducted for n E {1, 2, 4, 6}. 



III. NOTATION 

The 2D curved surface is embedded in Euchdean 3-space. If cr is a two-dimensional 
coordinate of a point P on the surface then -R(cr) is the position vector of P in the em- 
bedding space relative to some origin. Coordinate-dependent unit basis vectors Sa, where 
a takes the values 1 and 2 to label the coordinates, are defined in the directions of daR- A 
coordinate-invariant metric tensor is defined on the Riemannian 2-space, with components 
Qab = daR ■ dbR, from which the area of the closed surface is calculated: 



A 



I d'cT^ (3.1) 



where g = det g. Indices are raised and lowered by the metric and its inverse, in the 
standard manner. The extrinsic curvature tensor has components Kab = N ■ dadbR, where 
A?" is a unit vector normal to the surface. (N.B. Dot products are evaluated in M^.) Its trace, 
K^, is the 'total curvature'. Readers may be familiar with the mean curvature which is half 
of this quantity. The intrinsic or Gaussian curvature is given hj K = det I^. 

As an example, the above quantities are calculated for a sphere of radius R, using 
spherical polar coordinates {6, 0). The metric tensor may be represented by a matrix: 

( Q^^ Q^2\ J I ^ 
= R 



yO sin^^ J 



911 912 

y 921 922 j 

with determinant g = R'^sm'^9. Hence the area element for integration is ^/gdPcr = 
sin 6 d6 d(j), which is just the Jacobian for spherical polars. Any continuous, doubly dif- 
ferentiable surface has two principal radii of curvature, Ri and R2 at each point. These are 
the radii of curvature of the geodesies passing through that point which have respectively 
the greatest and least curvature as measured in the Euclidean embedding space. It can be 
shown [H] that these two geodesies are always perpendicular. If the particular coordinate 
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basis used corresponds to these two perpendicular directions, then the curvature tensor will 
be diagonal, and given by 



K = 



Kl Kl 



Kl Kl 



Of course, on a sphere, curvature is constant. Therefore K\ — where the sign convention 
denotes the direction of the normal. However, it should be noted that, without one index 
raised, the components are not constant: 



K21 K22 



-R 



^ ^ 



sin^^ 



Clearly, on a sphere, mean curvature is — and Gaussian curvature is 

Recall that the intrinsic n-atic order in the system is orientational, not positional, and 
that it is locally invariant under rotations through As this orientational order is intrinsic 
to a 2-dimensional space, in can be represented by a complex order parameter in the following 
way. Let 6 be the local angle between ei and the orientation of ordering. For vector order, 
this orientation is easily defined. For instance, in a smectic-C membrane, it is the direction of 
a molecule's aliphatic tail projected onto the local tangent plane of the membrane. However, 
in a nematic liquid crystal, may take one of two values for each molecule, depending on 
which end of the molecule one chooses to measure from. For our purposes, it does not matter 
which of these two values is assigned to at each point. For bond-angle orientational order, 
eg. hexatic, there is an n-fold ambiguity in the definition of but again, it does not matter 
which of the n orientations is chosen. It is conventional in this case to define at a given 
molecule as the angle between the basis vector ei and the line joining that molecule to its 
nearest neighbour. Having defined 0, the order parameter is given by 



■{/^(cr) = (exp(m0(cr))) 



(3.2) 



where <> denotes a thermal average. Note that ip has the appropriate rotational invariance, 
since a rotation of the orientation of ordering through an angle x within the surface is rep- 
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resented by a rotation of through nx in the Argand plane. Hence rotating the orientation 
of ordering through — corresponds an identity transformation in the complex plane. 



IV. THE MODEL 

The model used is identical to that in |^. There is zero spontaneous curvature as the 
bilayers considered are symmetric, and no constraint on bilayer area difference. The vesicle's 
volume is also unconstrained since small solvent molecules may diffuse freely in and out, so 
long as the solvent is very pure; containing no free ions, which would lead to osmotic pressure. 
The membrane's area A is fixed by the number of constituent molecules. The free energy 
functional 

F[^{a), R{a)]/T = J d'a^ (r|^|2 + \u\,^\^ + CD^^l,Dir + HK? + ^gK) (4.1) 

closely resembles the Ginzburg-Landau free energy for classical superconductors. For a 
discussion of how the thermodynamic coefficients relate to real physical and chemical quan- 
tities, see for instance 0. As defined in Eq. (|3.2|) , the phase of ip depends on the choice 



of basis vectors. Hence the derivatives Da in Eq.(O) must be covariant in order for the 
quantity D^-ipDl^ip* to be physical. Let the basis vectors be rotated through an angle — x(cr). 
Consequently, 

9ab — > g'ab = ^a^b9cd 

where A^(cr) are the components of the rotation matrix. Given that Da = da ~ inAa, 
demanding that D'^ipD^tp* is invariant under the gauge transformation implies that the 
components of A transform like 

Aa ^ A'a = A^iA, + d,a) . (4.2) 

So A undergoes the usual transformation, da<y, of a vector potential in a gauge-invariant 
derivative. But, in addition there is a coordinate rotation, A^ since, in this complex number 



representation, a rotation of the Argand plane is always accompanied by a rotation of the 
basis vectors. The 'spin connection', Aa = ei ■ dae2 has the appropriate transformation 
property Eg. ( [4.2| ). On the sphere, Ag = and A^ = — cos 6. 



V. THE CURVATURE ENERGY TERMS 

Above the continuous transition to the disordered phase, ip vanishes and only the last 
two terms in the free energy Eq.( [4.1| ) survive. These constitute the Helfrich Hamiltonian, 
describing the elastic energy of a constant area membrane. 

The intrinsic curvature energy, being proportional to the integral of Gaussian curvature, 
is a topological invariant, according to the Gauss-Bonnet formula 

jKdA = 2nx (5.1) 

where x is the Euler number of the closed surface [Q. For an orientable surface, x = 2(1 — G), 
where G is the surface's genus; ie. how many 'handles' it has. Hence the intrinsic curvature 
energy is zero for a torus and Attkg for a sphere. 

The total curvature energy, ^ J{K^)'^dA, is a con/orma/invariant [^j. ie. it is invariant 
under shape changes of the form 

R^R'= , (5.2) 

where /3 is a constant vector which parameterizes the conformal transformation . A sphere 
has total curvature energy Sttk and its shape is not changed by such a transformation. This 
is not true for tori. An axisymmetric, circular cross section torus may be deformed by a 
conformal transformation. As a consequence of the circular cross-section, a component of f3 
parallel to the symmetry axis serves only to change the torus's size. Hence there remains 
only a one-parameter family of non-trivial conformal transformations. With /3 lying in the 
symmetry plane, the surface is deformed into a non-axisymmetric torus such as the one in 
Fig.(0), and as (3 increases in magnitude, the asymmetry increases until, in the limit, the 
torus becomes a perfect sphere with an infinitesimal handle. 
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FIGURES 




FIG. 1. The surface produced by a conformal transformation of an axisymmetric, circular 
cross-section torus. 

The handle, although infinitesimal, must contain an intrinsic curvature energy of —Atckg 
and finite total curvature energy. A circular, axisymmetric torus of aspect ratio (ratio of 
radii of generating circles) a has total curvature energy 

^/,.C)=.-4^^, ,5.3) 

which is a minimum at a = , the Clifford torus, and infinite at a = 1, where the torus 
becomes self-intersecting. Willmore conjectured that Clifford's torus and its conformal 
transformations have the absolute minimum curvature energy of any surface of toroidal 
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topology; a conjecture that is difficult to prove, as the most general toroidal surface can- 
not easily be parameterized: tori can even be knotted. Accordingly, this paper will not 
explore the whole of phase space. The discussion will be restricted to axisymmetric, circular 
cross-section tori and their conformal transformations, with the assumption that these are 
close to the ground state shapes. To be consistent, the ground-state genus-zero shape is 
approximated by a sphere: the deviations from sphericity calculated in [§] are not included. 

VI. THE LOWEST LANDAU LEVEL APPROXIMATION 

The field gradient energy term in Eq.( |4.1| ) may be re-formulated. Using integration by 
parts, it is simple to show that, on a closed surface, 

J d^a^g'^'DaijDir = -J d'arD^i^g^'D,^). (6.1) 

Let 0p be eigenfunctions satisfying the Hermitian equation 

— = A0P (6.2) 

V9 

and the orthonormalization condition 

J d^a^4>l(pq = ^^pq- (6.3) 

These are the Landau levels. ip may be expanded in this complete set of functions. However, 
close to the mean-field transition, it is a good approximation to use only the Landau levels 
with the lowest eigenvalue, Aq. As is common in such problems, there is a degenerate set of 
lowest Landau levels. 

A scales inversely with area. So let us define a scale-invariant eigenvalue A, by 

An 

A = -A. (6.4) 

The free energy becomes 

F/T = y" I + 1^1^ + i'^'"^} ^ curvature energy 
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where 



(6.5) 



and Qp are complex constants. Note that the mean field transition has been lowered from 
r = to r = Tc where = — izE^Aa, Qn an intrinsically fiat surface, Aq = but it is finite 
and positive-definite on a surface with finite Gaussian curvature. So ordering is frustrated 
by intrinsic curvature. 

VII. SOLUTION OF THE EIGENFUNCTION EQUATION 



Eq. (|6.2| ) may be re- written as the differential equation 



..ab 



, db — 2inAij 
~ V9 



)iV99 



da4>p + 



-inD 



V9 



This is still general to any surface. 



A. Sphere 



Air 

ijp = -^A0p- (7-1) 



On a sphere, in spherical polar coordinates {9,(f)), Eq. (|7.1| ) becomes 



' g2 g2 g g N 

7— + cosec^6'7— ■ + cot 9— + 2incosec9 cot 9— cot^ 9 

09^ 09 00 ; 



which is solved by Aq = n and 0p of the form 



bp = -^<Pp (7.2) 



bp oc sm 



n+j / e 



cos 



(7.3) 



where p takes integer values from —n to n. These functions form an alternative basis set 
to that used in and are linear combinations of that set. They are given in |^ in a 
different gauge. The functions have 2n zeros, which are topological defects of index ^. The 
{2n + l)-fold degeneracy reflects the freedom of the defects to lie anywhere on the surface 
at no energy cost to order ip'^. 



11 



B. Axisymmetric Torus 



On tori, 0p and A depend on a and j3. Consider first the axisymmetric case {j3 = 0). Let 
us define coordinates {ri,(f)). The azimuthal angle is 0, varying around the large generating 
circle and subtended at the centre. The angle t] varies around the smaller generating circle 
and is measured from the symmetry plane, as shown in Fig.(^. 




FIG. 2. Ttie angular coordinates {r], (p) defined on an axisymmetric, circular cross-section torus. 
For clarity, only half of the torus is shown. 



Let the small generating circle have radius Ro- In this coordinate system, the following 
quantities are found: 



dab — Rq 







(a — cosriY ^ 



K. 



b ■ 



1 (-1 



a p 



g COS rj 
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— sm Tj ^ 

The Landau level equation Eq.( [7.1|) becomes 



dndn + 



+ 



'a — cosriY a — COST] 



-dr,+ 



2in sin rj 

[a — cos 



n 



sm ?7 



a — COST] 



AttX ^ 



(7.4) 



which is manifestly invariant under rj — >■ —rj with complex conjugation, reflecting the up- 
down mirror symmetry of the torus. Note that changing variable from rj to 9 = [t] — ^) and 
setting a to zero re-creates the sphere equation (Eq.( [7.2|) ). 

Eq. (|7.4| ) has separable variables and is solved by functions of the form 



MvA) = fpiv)e 



ip<f> 



(7.5) 



for integer p, where 



+ 



smrj 



d 



p + n sin f] 



drf (a — cos rj) dr] \ a — cos r/ 



47rA 



(7.6) 



and fp is real. It is immediately apparent from the form of Eq. (|7.6|) that 

fpiv) = f-pi-v) 



(7.7) 



and consequently 0p and 0_p are degenerate pairs of states. 

The ground-state solutions of Eq.( [7.6|) were found numerically for various values of p 
and a, using Rayleigh-Ritz variational method. Trial functions were synthesized from a 
truncated Fourier series, 

^max 

/p = i^^ + ^^"^ ' 

which has the appropriate periodic boundary conditions. The Fourier coefficients were used 
as variational parameters to minimize the integral 



J^\a-cosr])f;-' 



d'^ sin T) d I p + n sin rj 

dr]"^ {a — cos rj) drj \ a — cos r] ^ 



(7.8) 
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subject to Eq.( |6.3|) . The minimal value of Eq.( |7.8| ) is an upper bound on Aq , and the 
corresponding function /p""*' is approximately fp . Using a Fourier series truncated after 
frequency u,-,-,^^ creates a {21/,-^^^ + l)-dimensional parameter space to explore in the numerical 
minimization. So CPU time increases very rapidly with u^^^. A second numerical method 
was also used, which involved numerically 'integrating' Eq.( [7.6D forward in 77 from rjQ to 
riQ + 2tt and minimizing the resultant discontinuities in fp and its derivative, with respect to 
the initial conditions and A. This method is more accurate than the Rayleigh-Ritz method 
but more difficult to use since, depending on the starting conditions, it may find a non- 
ground-state solution. However, it was useful in establishing the accuracy of the upper 
bounds. Using u^^^ = 7, Aq was found to an accuracy of 0.05% in the worst case. 

In Fig.(|]), A(a) is plotted for < p < n, with data points at increments of 0.5 in a 
and cubic spline interpolation. The states are labelled Up. One might have guessed that uq 
would be the ground state in each case. And indeed, for large aspect ratios, when the torus 
is almost cylinder- like, the p = states do have the lowest energy. These states are even 
functions of rj, being larger at rj = it than aX rj = where there is more absolute intrinsic 
curvature, and do not vary with 0. ie. The field does not twist with respect to the coordinate 
lines as (p increases. However, for smaller a, states with non-zero p are energetically favoured. 
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FIG. 3. Graphs of A against a for axisymmetric, circular cross-section tori. The states are 



labelled Up. 



Before explaining the reason for this, it is necessary to describe these states. Fig.(^ 
depicts a selection of states for = 2 on a Clifford torus. The nematic-order field is 
represented by sets of double-ended lines directed along the orientations of order in the 
surface of the torus. The length of each line is proportional to {ipl at that point. Note however 
that the three-dimensional views of the toroidal surface cause the lines to be foreshortened 
at obhque angles to the hne of sight. 
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FIG. 4. Various Landau levels of the n = 2 (nematic) order-parameter field on a Clifford torus. 
States are labelled Hp and shown in plan view and side view and 77 dependence is given graphically. 
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Two views are shown of each field on the toroidal surface, as well as a graph of fit]). 
Note that positive p states have a large field on the bottom of the torus and small on the top, 
and vice versa for negative p. This asymmetry becomes more pronounced with increasing 
p, with increasing n, and with decreasing a. As stated earlier, ordering is frustrated by 
intrinsic curvature. And the net amount of intrinsic curvature on the torus is constant. But 
as a decreases, the curvature is redistributed. It becomes more negative around the hole and 
more positive around the outside of the torus, with the top and bottom remaining fairly flat. 
So this is where the field prefers to reside and be parallel. Although in Eq.( [r.5|) non-zero 
p states appear to twist with respect to the coordinate basis as (p increases, the basis itself 
twists on the top and bottom of the torus as the hole is circled. The two twists cancel. 
Note that the diagram of the 2_2 state depicts a field which is as parallel as possible on the 
top 'face' of the torus. Of course, the field on the inner and outer 'rims' still contributes 
to the free energy so, as a increases, the ground-state value of p decreases by a series of 
'transitions'. This is clear in Fig.(^. These should not be regarded as true phase transitions, 
since the lowest Landau level approximation breaks down completely at these points. As a 
crossover between eigenvalues is approached, the true state is a mixture of the two states. 

It should be emphasized that the tori studied here are true, curved tori, embedded in 
3-space. They should not be confused with a fiat torus, which is a plane with opposite edges 
identified, for which the above energetics would not apply. The two have the same topology, 
but different local geometries. 



When /3 7^ 0, the Landau level equation (Eq. ( |7.1|) ) becomes much more complicated and 
its variables are no longer separable. However, Rayleigh-Ritz variational method is again 
possible, as 



C. Non-Axisymmetric Torus 




(7.9) 
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but the trial wavefunctions are now complex, periodic functions of two variables. Thus, syn- 
thesizing them with a truncated, two-dimensional, complex Fourier series up to frequency 
^'max produces an [Suf^^^ + Su,^^^^ + l)-dimensional space of variational parameters to be nu- 
merically explored. Hence ground-state functions of the coordinates with high frequency 
components cannot viably be found. This is unfortunate since conformally transforming an 
axisymmetric torus with the coordinate system described above, leaves the coordinates very 
unevenly distributed over the surface, with most of the interval from to 2tt in rj and 
covering only the small 'handle' part of the torus and only a small range covering most of 
the area. This problem is partly solved by a coordinate transformation [rj, 0) —>■ (p, a) which 
preserves the metric's diagonalization, where 

[1 - /3(a - 1)] tan ^ = [1 
[1 + /?(« + 1)] tan I = [1 

although this still leaves some high frequency components in the distribution of the coordi- 
nates on the handle. Using these coordinates, Eq.( |7.9| ) is very complicated but two special 
cases are already known: 

• When j3 = 0, the torus is axisymmetric and Aq was found by the method in section 

vitbi . 

• When P = , where Ro is the radius of the small generating circle, the torus 
becomes a sphere with an infinitesimal handle. Hence Xq = n and the ground-state 
function is that of a sphere with all 2n zeros at the handle. 

Unfortunately, the problem is too computer-intensive to produce any useful results for inter- 
mediate values of (3. Using u^^^ = 3 and performing the minimization in a 98-dimensional 
parameter space took a considerable amount of CPU time per data point. At this level, no 
states of intermediate (3 were found to have a lower eigenvalue than the known axisymmetric 
or spherical states. Some local minima were found but these are very likely to be artifacts 
of the varying reliability of the upper bound, which was high by around 10% for /3 = and 
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— j3{a + 1)] tan 



V 



f3{a + 1)] tan — 



60% in the spherical hmit. More research on this problem is planned. But for the moment, 
it is necessary to work with the ansatz that: 

• A spherical vesicle exists when its total free energy is lower than that of an axisym- 
metric torus. 

• An axisymmetric torus exists when its total free energy is lower than that of a sphere 
and its field energy alone is lower than that of a sphere. 

• A non-axisymmetric torus exists otherwise. 

VIII. THE ip'^ TERM 

A. Axisymmetric Torus 

It was noted earlier that (f)p and 0_p are degenerate to quadratic order in ip. So, from 
Eq.( |6.5D , is a linear combination of (pp and where p has the ground-state value for 
the given aspect ratio. The relative proportions of the two states are determined by the ip'^ 
term in Eq.dOp. Let 

= a^(pp + a_(p^p 

be the function which minimizes / |\E'|^(i^ on a surface of area An subject to |a+p + |a_p = 1, 
and let J be the minimal value of that integral. Then Eq. ( ^4.1|) is minimized hj ip = a'^ in 
the regime of the lowest Landau level approximation, where a is real. The gauge invariance 
and axial symmetry of the problem allow us to choose the phases of both a+ and a_ . Hence 
is determined by a minimization with respect to just one parameter. 
For some values of a and p, it is found that \E' = + ct-)- So the 'top' and 'bottom' 

states co-exist and, although each is fairly parallel on its own side of the torus, the two 
do not mesh smoothly where they meet. Consequently, there are 2p defects each of index 
- around the torus's outer 'rim' and 2p of index — round the inside of the hole: a total 
of 4p defects. Hence as a varies and p jumps between integer values as described above, 
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vortex-antivortex quartets are formed or destroyed. Let these symmetrical states be labelled 
rip. The 4f state for a = \/2 is represented in Fig.(||). 



/ '^^ ^ \ f f ^ 



ixx X X X X X X X X ;:■ j UI 



^^"^ x">^ >^ % ^ V \ 



■x~ •><- ^ 
^ X. X X X XXXx 

X X X X X \\\ 




FIG. 5. 'Tetratic' {n = 4) order-parameter field in the 41 configuration. The field is represented 
by crosses of four lines, directed along the orientations of order in the surface of the Clifford torus. 
The length of each line is proportional to |^| at that point although the three-dimensional view 
of the toroidal surface causes the lines to be foreshortened at oblique angles to the line of sight. 
Notice that 3 of the 6 defects of index ^ are visible in the side view (where lip] ^0), and that the 
field rotates through = | of a turn around the hole. 

For other values of a and p, where the overlap between 0p and 0_p is large and hence 
the defects in the field become expensive, the up-down symmetry is spontaneously broken, 
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so that \1/ = 0p or \1/ = 0_p. Let these broken-symmetry states be labelled rip. As a varies, 
there are transitions between S and B states, which, unlike the transitions betwen states 
of differing p described in section |V11 B| , may be regarded as true phase transitions within 



mean-field theory, although they may be blurred by the inclusion of fluctuations. As a rule of 
thumb, the ip'^ term chooses S or B on the basis of which covers the surface more uniformly. 

A table, giving the values of a for all transitions between ground-state field configurations 
on the torus, is produced in the appendix. 

The free energy is next minimized with respect to a, giving the overall magnitude of the 
order parameter field. Below the continuous transition, the free energy is a minimum when 

_ / _l_ 47rCAo \ 

a' = ^ / K (8.1) 
uJ 



B. Sphere 



It was explained in section |VII A| that the spherical Landau levels are degenerate with 
respect to the positions of the vortices. The ip"^ term lifts this degeneracy and makes the 
vortices repel. Hence, following Park et al. [^, the ground state configuration is taken to 
be that of highest symmetry. So, when n = 1 the two zeros are antipodal. When n = 2, 
the zeros lie at the corners of a regular tetrahedron and, when n = 6, at the corners of an 
icosahedron. The values of J(sph){i^) were found using these configurations. For n = 4, the 
defects lie at the corners of a cube whose top face has been rotated through 45° relative to 
the bottom face. The distance between these two faces is then varied to minimize J. The 
optimal separation is found to be very close to the length of a side of the original cube. 

IX. THE PHASE DIAGRAMS 

The minimized free energy of the sphere is now completely determined. Expressing 
it as a fraction of the bending rigidity gives an equation relating just four dimensionless 
combinations of thermodynamic parameters: 

21 



Tk 



An (2 + 



1^ 



- 2n 



Auk 



Ar 

C 



+ n 



Jisph){n) 



for r < 
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A ' 

otherwise. (9.1) 



On the torus, the free energy remains a function of a only: 

27r 



F(tor) 



2'n'a 



2^,2 



Auk 



J(^,,){a,n) 



for r < 
otherwise 



-47rCAo(a) 
A 



(9.2) 



which relates just three dimensionless combinations of thermodynamic parameters. When 
Eq.( p.2|) is minimized with respect to a, some of the physical richness of the problem is lost, 
as many of the interesting phases and all the broken symmetry states occur on tori of aspect 
ratios which are not the energetic optimum. However, other interesting features arise, as 
will be seen in the phase diagrams in Figs.(|^)-(^), which result from a phase space restricted 
to only axisymmetric tori, and show contours of constant a as well as lines of first order and 
continuous phase transitions. The phases are labelled with the notation introduced earlier. 
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FIG. 6. Phase diagram for a phase space restricted to axisymmetric, circular cross-section tori 
of constant area A and variable aspect ratio a. Vesicles in the low-temperature phase have intrinsic 



vector (n = 1) order. See Eq.(4.1) for definitions of the thermodynamic parameters. 
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FIG. 8. Phase diagram of axisymmetric, circular cross-section tori with intrinsic 'tetratic' 
= 4) order in the low temperature phases. 
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J? r 

FIG. 9. Phase diagram of axisymmetric, circular cross-section tori with intrinsic hexatic (n = 6) 
order in the low temperature phases. 

The diagrams are plotted in the plane of the two independent combinations of ther- 
modynamic parameters in Eq. (|9.2|) . In each case, to the far right of the diagram is the 
high-temperature phase, for which ip = 0. In this phase, only the curvature energy is non- 
vanishing so the aspect ratio is \/2; the Clifford torus. Note that the continuous transition, 
at which ordering arises, occurs not at r = 0, but is suppressed by the manifold's curvature 
to r = Tc, as defined in section |V1| for the Clifford torus. Note that depends on aspect 
ratio, aspect ratio depends on ip and depends on Vc- This non-linearity is the cause of 
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the first-order phase transitions. Consider for instance the n = 1 diagram (Fig.(^). A 
vesicle with small k is easily deformed by the 1q -state ordering to large aspect ratios at 
low temperatures. At large aspect ratios, Vc is close to zero. As the temperature rises, ip 
diminishes, and can no longer hold the aspect ratio so far from the Clifford torus. As the 
aspect ratio falls, ordering becomes more difficult. There comes a temperature at which 
the torus jumps down a first order phase transition to the Clifford torus, where Tc < r. At 
higher values of k, the transition is continuous. The n = 1 phase diagram has a tri-critical 
point at ^ = —13.68 and ~ 0.98. For higher values of n, the first order transition 
persists for all values of k, and meets the continuous transition at a multi-critical point. 
The continuous transition here is between the high-temperature phase and the n^_i phase, 
for which a < a/2 since the gradients in Fig.(^ are positive close to the Clifford torus. 

The free energy (Eq.( |4.1|) ) has now been minimized with respect to the order parameter 
field's configuration and magnitude and with respect to the toroidal aspect ratio. It remains 
only to minimize with respect to topology, according to the prescription given in section 
|V11 C| . Figs, (pup -(]TB|) show the resulting phase diagrams, which include spherical, and ax- 
isymmetric and non-axisymmetric toroidal phases of vesicles. No a contours are shown in 
these figures. The positions of the transitions vary with and transition lines are shown 
for various values of this quantity. Figs. ([ToD -([l3|) should therefore each be regarded as several 
phase diagrams superimposed. 
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FIG. 10. Phase diagram for constant area vesicles of genus zero or one, with intrinsic vector 
(n = 1) order. The position of the transition hne between genera depends on the ratio of intrinsic 
and total bending rigidities In the un-shaded region, genus one vesicles are axisymmetric, 
circular cross-section tori, whose aspect ratio is given in Fig(P). Genus one vesicles that exist in 
the shaded region are non-axisymmetric. The area of the diagram occupied by spherical vesicles 
decreases with increasing 
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FIG. 12. Phase diagram for constant area vesicles of genus zero or one, with intrinsic tetratic 
= 4) order. 
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FIG. 13. Phase diagram for constant area vesicles of genus zero or one, with intrinsic hexatic 
(n = 6) order. 

The phases do not persist when this larger phase space is explored, except for 6f , 
for which Aq is less than the sphere value of 6, when a is close to ^/2. The diagrams are each 
divided into shaded and unshaded regions by a heavy line. In the shaded regions, vesicles 
are either spheres or non-axisymmetric tori while, in unshaded regions, either spheres or 
axisymmetric tori occur. 

Consider first the n — 1, n — 2 and n — A phase diagrams. A transition line, whose 
position depends on ^ separates tori from spheres. Note that for ^ > tt — 2, the transition 
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line is closed, with spheres on the inside and tori of one sort or another on the outside, while 
for ^ < vr — 2, the line is open, with spheres on the high-temperature side and axisymmetric 
tori on the low-temperature side. This magical value of vr — 2 is the critical value of ^ at 
which spheres would become tori if the model described curvature energy alone and did not 
include terms in tp. Note that in these three phase diagrams the separatrix (the ^ = tt — 2 
transition) has a vertical limb at the value of ^ at which ordering first forms on the sphere, 
le. at r = r=p'^'='°. 

The n = 6 diagram in Fig. ([13|) is a little more complicated since there are two distinct 
phases of ordering on axisymmetric tori. Again, the diagram is divided into shaded and 
unshaded regions. The shaded region is now split into two parts. The part on the right 
contains spheres if ^ < vr — 2, and non-axisymmetric tori otherwise. The part on the left 
contains spheres on the high side of the transition line (whose position depends, of 
course, on ^) and non-axisymmetric tori on the other side. The un-shaded region, which 
contains axisymmetric tori and spheres, is also divided into two parts. The part in the top 
left of the diagram behaves as described above for the other phase diagrams. The other 
part may (depending on ^) contain a closed transition line, which, unlike the other closed 
transition lines we have met, has spheres on the outside and axisymmetric tori, of small 
aspect ratio and 6f field configuration, on the inside. 

The four phase diagram of Figs. (|T0D -(p!3D may seem complicated but they are easily 
understood if one remembers that: a) they are several diagrams superimposed, and b) 
the regions in which spherical vesicles exist become smaller as increases, since spheres 
have more intrinsic curvature than tori. Having established from these diagrams that an 
axisymmetric torus exists in a certain region, its aspect ratio may be found from Figs.(^-(^. 

It should be mentioned that the curves defining the edges of the shaded regions in 
Figs.(|T0|)-(|T3D have been slightly smoothed by hand to remove pixcellation which arose from 
numerical rounding errors. 
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X. VALIDITY OF THE APPROXIMATIONS 



Some order-of-magnitude calculations are now presented to give an indication of the 
regimes of validity of the various approximations used. 

A. Lowest Landau Level Approximation 

The Landau levels form a complete set of orthogonal functions in which to expand 
the order parameter; ip = J2p(ip(pp- Each eigenfunction 0p has an eigenvalue Ap defined 
by Eqs .( 16.21) and (|6.4|) . The lowest Landau levels alone may be used if <^ Og where 
Ap 7^ Ag = Aq. As an indicator of the region in which this is valid, let us calculate < a*ap > 
for Gaussian fluctuations above the continuous transition, using the quadratic Hamiltonian 

H = ^(r-7Ap)6;6p + 0(6^) (10.1) 

p 

where 7 = — and bp = \f^ap- The partition function is 

2 = exp^(7Ap - r)blbp Hdb^dbp. 
p p 

Hence 

So for the lowest landau levels, < b*bp >o= and for the next lowest level, < b*bp >i = 
Hence the ratio of correlators is 



r— 7A 

7^ 



< a*pap >i _ < b*bp >i _ x 

< a*ap >o < b*bp >o x + 47r(Ai - Aq) 



where x = ^^ is the x-coordinate in the phase diagrams, measured with respect to the 
transition. We see that TZ is small, and hence the lowest Landau level approximation is 
good, when |a;|<47r(Ai — Aq), the right hand side of which may be estimated from Fig.(^. 
For the Clifford torus, (Ai — Aq) ~ 0.5 so the approximation is quantitatively accurate only 
within ~3 of the transition, along the x-axis. However, Figs.(^-(|l^) cover a much larger 
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region than this, in order to clarify the quahtative features of the diagrams, for which the 
calculation at this level of approximation is a useful guide. It should be noted that, at 
'transitions' from one lowest Landau level configuration to another, (Ai — Aq) 0, so the 
approximation must break down here, as noted in section |V11B| . Also, in the large-a phases 
below the first-order transitions, is higher, so the system is far from the lowest Landau 
level regime. These phases are included in the diagrams in order to show the first-order 
transitions. 



B. Mean Field Theory 

For this order-of-magnitude analysis, let us treat two parts of the free energy separately: 
first we will demand that fluctuations in the Helfrich Hamiltonian may be neglected in the 
absence of orientational order, then, within this fixed, non-fiuctuating environment, it will 
be required that the terms in ip satisfy the Ginzburg criterion for the validity of mean-field 
theory. 

In a simplistic calculation, good over only a short range, the Monge gauge is used to 
parameterize a locally horizontal surface, defining its z-coordinate as a function of x and y. 
Keeping terms of second order in z, and expanding z in plane-wave modes: z ~ Y,k cte^^'^ 
gives 

K / {KfdA ^ K [ {V'^zfdxdy = nA'^Y.Wkl^k^- 

Note that the Gaussian curvature energy does not influence fluctuations as it is a topological 
invariant. Hence 

1 



nA^k^ 

so, for fluctuations to be negligible on the scale of the system size requires k 3> ^ . This 
is intuitively obvious; that a large bending rigidity is required to combat fluctuations. In 
fact, this simple calculation on a flatish surface is over-stringent, since intrinsic curvature, 
as exists in closed surfaces, reduces the size of fluctuations [0. 
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The Ginzburg criterion is now computed for the Ginzburg-Landau part of the free en- 
ergy. This is a comparison between the size of the discontinuity in specific heat at the phase 
transition, calculated with mean-field theory, and the size of the divergence in specific heat 
due to Gaussian fluctuations. Close to the transition, r goes approximately linearly with 
temperature T. Hence (r — r^) = p(T — Tc), where is the shifted mean-field transition tem- 
perature, and p is some constant. From Eq.( |9.1|) it follows that the mean-field discontinuity 
in specific heat is 

For simplicity, the specific heat of fluctuations is calculated for the sphere, above the phase 
transition. Using the quadratic form Eq. ( |10.1| ), and only the (2n+ 1) lowest Landau levels, 
gives a free energy density / = ^^"^^^"^ (In p(T — T^) — Iutt). Hence the singular part or the 
specific heat is 

_ (2n + l)p 

^sing. 



giving a critical exponent a = 2 since we are in a zero dimensional fc-space. The ratio of 
these two specific heats is 

C If 

^sing. / S 



AC U-f 



where 



c 



/47r(2n + 1)mJ 



/ ApT, 

and mean field theory is good for |r— rc|>C • (N.B. This condition is equivalent to demanding 
that tree graphs dominate over one- loop graphs.) So the threshold, in the phase diagrams, 
between the mean-field regime and the fluctuation-dominated regime is given by the curve 
x^y = (2n + where x = -^^^'^""^ and y = having set p ~ T^^ and J ~ At 



a height yo, this region has width A = y ^ ^yo • '^^^^ should be much smaller than the 



region of interest ie. that in which the lowest Landau level approximation is valid. So A ^ 3 
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and, by inspection, in all the phase diagrams, we are interested in ?/o ~ ji- Consequently the 
condition for a mean-field treatment of the order parameter to be valid in this investigation 

2 

is just K ^ which is always consistent with the validity of a mean-field treatment of the 
membrane shape, calculated above. 



C. The Ginzburg-Landau Model 

The Ginzburg-Landau free energy functional is and expansion in ip, truncated at order 
?/'^. This is valid when vAI'j/'P <^ 1. In the mean-field regime, this is always the case above 
the continuous transition. Below the transition, ijj ~ from Eq. (|8.1|) . So the inequality 



uJv 47r 

is marginal when y ~ ^^^^ where x and y are given above. This defines a region in the phase 
diagrams of width A = -^j^ at a height y = ^, which is in the region of interest. Again, 
let's demand that A >> 3 to be consistend with the region of validity of the lowest Landau 
level approximation. Consequently, it is required that § ^ 

Dropping constants of order unity, we find that mean-field theory, the Ginzburg-Landau 
model and the lowest Landau level approximation are all valid in a region of half- width ~ 3 
in the phase diagrams, given that 

— > K > 1. 
A 



D. The Shapes 

Close to the ordering transition, ifj is very small, so the order cannot deform the torus 
far from Clifford's aspect ratio \/2, or the genus-zero shape far from a sphere, so the vesicle 
shapes used in this theoretical investigation are undoubtedly good. As Park et al. discovered 
in 0, wants to live in a fiat space, and can achieve a fiattening of the membrane via the 
gauge-like coupling in Eq.( |4.1| ). This it will do where it is large, at the expense of an increased 
curvature where ip is small eg. at the defects. This pay-off is unavoidable, since the integral 
of intrinsic curvature is topologically fixed. Hence spherical vesicles are slightly deformed 
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by the in-plane order, into rounded polyhedra, whose vertices are at the order parameter's 
vortices. This is the effect which has been ignored in this paper. The same process will also 
occur in tori, especially in those parts of the phase diagrams where ordering is predicted to 
deform the tori to large aspect ratios. In general, one might expect that the larger the aspect 
ratio, the worse the approximations of circular cross-section and axisymmetry. But nature is 
on our side here. Since, in these large-aspect-ratio tori, if) is always in the Uq configuration, 
which has no vortices, and becomes very uniform for large aspect ratios. So, although this 
approximation results in some quantitative inaccuracies, they are likely to be small. Also, 
they must change the free energy of the torus and the sphere, if not by the same amount, 
then at least in the same direction. So the phase diagrams may not have exactly the right 
numbers, but they do have more-or-less the right shape. 

XI. CONCLUSION 

It is clear that orientational order in a toroidal geometry is a subject rich in physical 
diversity, and deserving of further study in condensed matter theory. Also, introducing 
orientational order into the surface is an interesting way of stabilizing a torus against the 
massless Goldstone mode of conformal transformations, without resorting to spontaneous 
curvature or a constraint of constant volume or bilayer area difference. As far as experimental 
verification is concerned, the principal results of this paper are embodied in the final four 
phase diagrams. The approximations used in their derivation may result in quantitative 
inaccuracies, but the main qualitative features of these diagrams are predicted to be seen in 
experiment in the future. At present an experimental realization is somewhat difficult, due 
to the fragile nature of laboratory-produced vesicles. 
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APPENDIX A: TABLE OF MEAN-FIELD PHASE TRANSITIONS 



The axisymmetric toroidal aspect ratios at which the order parameter field, which min- 
imizes the free energy in the regime of the lowest Landau level approximation, changes its 
configuration between states of different p or different symmetry are given in the table be- 



low. See sections [Vll B| and |V111 A| for explanations of these transitions. Note that an n = 1 
(vector) field has no such transitions, as its ground state is always 1q . 



Transition 


Aspect Ratio 




2.31 
2.41 


4f ^ 4f 
4f - 4f 


3.34 
4.56 

4.98 


6f 6f 
6f ^ 6f 
6f - 6f 

6f 


3.15 
5.46 
6.63 
7.24 
7.83 



38 



REFERENCES 

[1] F.C.MacKintosh and T.C.Lubensky, Phys.Rev.Lett. 67 1169 (1991). 

[2] Jeongman Park, T.C.Lubensky and F.C.MacKintosh, Europhys.Lett. 20,3 279 (1992). 

[3] D.E.Rutherford, Vector Methods (Ohver & Boyd Ltd., 1959). 

[4] Jeffrey R. Weeks, The Shape of Space (Marcel Dekker, 1985). 

[5] B.Fourcade, J.Phys.II France 2 1705 (1992). 

[6] F.Jiihcher, U.Seifert and R.Lipowsky, J.Phys.II France 3 1681 (1993). 

[7] T.J.Willmore, Total Curvature in Riemannian Geometry, (Elhs Horwood Ltd., Chices- 
ter, 1982). 

[8] J.A.O'Neill and M.A.Moore, Phys.Rev.B 48 374 (1993). 

[9] F.David, Geometry and Field Theory of Random Surfaces and Membranes, in Statistical 
Mechanics of Membranes and Surfaces, ed. D.Nelson, T.Piran and S.Weinberg (World 
Scientific Publishing Co. Pte. Ltd., 1989) 

[10] R.de Vries, Private communication. 



39 



